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We present studies of the magnetic field distribution around the vortices in LuNi2B2C. Small- 
angle neutron scattering measurements of the vortex lattice (VL) in this material were extended to 
unprecedentedly large values of the scattering vector q, obtained both by using high magnetic fields to 
decrease the VL spacing and by using higher order reflections. A square VL, oriented with the nearest 
neighbor direction along the crystalline [110] direction, was observed up to the highest measured 
field. The first-order VL form factor, \F{qio)\, was found to decrease exponentially with increasing 
magnetic field. Measurements of the higher order form factors, \F{qhk)\; reveal a significant in-plane 
anisotropy and also allow for a real-space reconstruction of the VL field distribution. 

PACS numbers: 74.25.Qt, 74.25.0p, 74.70.Dd, 61.05.fg 



I. INTRODUCTION 

The magnetic field distribution due to the vortex lat- 
tice (VL) in type-II superconductors depends on the de- 
tailed nature of the superconducting state and on the 
properties of the host material. Examples of different 
field profiles are evident if one considers the results of 
calculations based on different theoretical models for the 
superconducting state.— —i^— Experimentally one often 
seeks to parameterize the field modulation in terms of two 
characteristic length scales: the penetration depth (A) 
and the coherence length (^) . Such an approach provides 
a simplified method of analyzing the results of small- 
angle neutron scattering (SANS)f^ muon spin rotation 
(/iSR)^ and nuclear magnetic resonance (NMR) measure- 
ments. However, in addition to the simplification, such 
an approach also requires the implicit acceptance of a 
particular theoretical model while, in many cases, vio- 
lating its premises by for example using the Ginzburg- 
Landau model to extract a field dependent penetration 
depth and coherence Icngthjii^ 

In this paper we will describe a more complete, model- 
independent analysis of SANS measurements of the VL 
in LuNi2B2C, extended significantly beyond the first or- 
der Bragg reflection, which is customarily the only one 
measured. Measurements of a large number of reflections 
allows for a real-space reconstruction of the VL magnetic 
field profile, which will be discussed in relation to the 
significant in-plane anisotropy of this material caused by 
the Fermi surfac o^'^^d^'^^ and the anisotropic pairing in 
the superconducting stateFi^ii^ii^ii^iiiii^ii^i^a To the best 
of our knowledge only very limited efforts have been un- 
dertaken in measuring higher order VL reflections, the 
most notable exception being the work on Sr2Ru04by 
Kealey et al^ 



II. EXPERIMENTAL DETAILS 

LuNi2B2C is a non-magnetic member of the rare-earth 
nickelborocarbide family of superconductors with a criti- 
cal temperature of Tc = 16.6 The single crystal used 
in the SANS experiment was grown by a high tempera- 
ture flux method)^ using isotopically enriched ^^B to re- 
duce neutron absorption, and subsequently annealed to 
improve quality and reduce vortex pinning The sample 
had a mass of ~ 1 g and a disc-like crypto-morphology 
with the c axis parallel to the thin direction. 

The experiment was performed at the Dll SANS in- 
strument at the Institut Laue-langevin. Incident neu- 
trons with wavelength A„ = 0.45 nm and wavelength 
spread of AA„/A„ = 10% were used, and the VL diffrac- 
tion pattern was collected by a position sensitive detec- 
tor. Measurements were performed at 2 K in horizontal 
magnetic field between 0.5 and 6 T, applied parallel to 
both the crystalline c axis and the incoming beam of neu- 
trons. Two different magnetic field-temperature histories 
were employed: Field cooled (EC) from a temperature 
above and zero field cooled (ZEC) followed by an in- 
crease of the magnetic field at 2 K. 



III. RESULTS 

Here we present SANS imaging of the VL in LuNi2B2C 
to an unprecedentedly high field of 5.5 T correspond- 
ing to 75% of the upper critical field, Hc2{'2 K) = 7.3 
rj 24,25,26 ^-|, fields a square VL was observed as shown 
in Eig. [TJa). An indexing of the VL Bragg reflections is 
shown in Eig. [ijb). The different reflections are located 
at a distance from the center of the detector which is pro- 
portional to their scattering vector, qhk = (/i^ + fc^)^/^ qo; 
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FIG. 1: (Color online) SANS diffraction pattern of the VL in LuNi2B2C at 0.5 T and 2 K following a field cooling procedure. 
The image (a) is a sum of measurements as the sample is rotated around the vertical axis in order to satisfy the Bragg condition 
for refiections in the center-right part of the detector. The data is smoothed and shown on a logarithmic scale. Measurements at 
6 T, where no scattering from the VL could be observed, were used for background subtraction. The azxes show the orientation 
of the crystalline axes. An indexing of the peaks is shown in (b), with the cross indicating the origin of reciprocal space. The 
apparent difference in intensity of e.g. the (12) and (21) reflections compared to (12) and (21) is due to different Lorentz 
factors, for which the detailed reflectivity analysis have been corrected. 



where qo = 2tt{B /(jjoY^^ and (/)o = h/2e = 2070 Tnm^ is 
the flux quantum. With increasing field, the VL Bragg 
peaks move out in reciprocal space and their intensities 
decrease, and as a consequence fewer peaks are visible. 
At 5 T and above only the {10}-reflections are observed. 
Measurements performed at 6 T were used for back- 
ground subtraction. While this is below Hc2 no scat- 
tering from the VL could be observed at this field. Fur- 
thermore, the detailed measurement in Fig. [4] show that 
at (jio(6 T) = 0.34 nm^^ the extrapolated (10) intensity 
is more than an order of magnitude smaller than that of 
any observed reflection at lower fields at the same q. 



flections broaden significantly in the azimuthal direction. 
While this could be due to the onset of recntrance of 
the square VL phase (transition back into a rhombic 
symmetry no splitting into two peaks was seen, and 
therefore the broadening may also simply be due to a dis- 
ordering of the vortex lattice. No indication of a 45° VL 
rotation was observed. Whereas this does not exclude 
such a transition at even higher fields it significantly re- 
duces the fraction of the HT-phase diagram where it 
can occur. This result therefore imposes significant con- 
straints on the model parameters used in the calculations 
described above. 



A. Vortex lattice symmetry and orientation 

At low applied magnetic fields the VL in LuNi2B2C 
undergoes a field-driven symmetry- and reorientation 
transitioni22i2^ This two-step transition arises due to the 
growing importance of the Fermi surface anisotropy cou- 
pled with non-local electrodynamics as the vortex den- 
sity increases)^ At higher fields, it has been proposed 
theoretically that thermal vortex fluctuations may lead 
to a re-entrance of the square VL phase.— A similar re- 
entrance was also predicted by Nakai et al. who consid- 
ered a case where competing superconducting gap and 
Fermi surface anisotropics both favor a square VL, but 
oriented at 45° with respect to one another.— In addition 
to the re-entrance of the square VL phase stabilized by 
the Fermi surface anisotropy, this model predicts a 45° 
rotated square VL phase at even higher fields due to the 
gap anisotropy. 

Experimentally wc found that above 5 T the VL re- 



B. Vortex lattice reflectivity and form factor 

Measuring the intensity of the VL reflections as the 
sample is rotated around the vertical axis to satisfy the 
Bragg condition, provides rocking curves as the ones 
shown in Fig [2l In addition to the strongest (10) re- 
flection, the figure shows the (32) rocking curve which 
was the highest order reflection visible at a field of 0.5 T. 
The intensity of these two reflections differ by a factor 
of 3000. The longer scattering vector for the (32) reflec- 
tion, = vTSqo, is evident by the larger rotation an- 
gle necessary to satisfy the Bragg condition. To obtain 
the VL reflectivity, the integrated intensity was deter- 
mined by fitting a Voigt function to each rocking curve 
and normalizing the area to the incident beam intensity. 
Compared to other functional forms (e.g. Gaussian or 
Lorcntzian) the Voigt was found to provide a significantly 
better fit to the data. The difference in the Lorentz fac- 
tor (angle between the scattering vector and the vertical 
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FIG. 2: (Color online) Rocking curves at 0.5 T and 2 K for the 
LuNi2B2C (10) and (32) VL reflections from Fig.[Tl Note the 
different intensity scales for the two reflections. Error bars for 
the (10) reflection are not shown since they are smaller than 
the size of the data points. The intensity at each angular 
setting are obtained by summing the detector counts at the 
position of the specific Bragg reflection. The curves are Voigt 
fits to the data. The shoulder seen for the (32) reflection is 
unrelated to the VL, as discussed in the text. 



rotation axis) for the two reflections gives rise to a dif- 
ference in the width of the two rocking curves, as they 
are cutting through the Ewald sphere at different inci- 
dent angles. The integrated intensity for all reflections 
is corrected for this effect. What ultimately limits how 
many VL Bragg peaks can be imaged is the vanishing in- 
tensity and imperfect background subtraction as seen for 
the (32)-rcflection in Fig. [21 Here significant background 
variation is clearly evident, even leading to an apparent 
shoulder on the VL rocking curve. While this could be 
interpreted as being due to a second VL domain, this is 
clearly not the case since a similar shoulder is not seen 
on the (10) rocking curve. 

With the strong VL peaks in LuNi2B2C, especially at 
low fields as shown in Fig. [1] it is necessary to consider 
whether multiple scattering is affecting the measured in- 
tensities. Four multiple scattering processes can affect 
the measured intensity of a given Bragg reflection: (i) Ex- 
tinction by scattering back into the incident beam, (ii) 
a diminished incident beam intensity due to scattering 
into other reflections, (iii) scattering into other reflections 
("aufhellung"), and (iv) scattering from other reflections 
( "umweganregung" ) The first 3 processes decrease the 
intensity of the specific rcfiection whereas the last one in- 
creases it. Furthermore for (ii) to (iv) to occur multiple 
VL reflections must satisfy the Bragg scattering condi- 
tion simultaneously. Since the effects of multiple scatter- 
ing depend on the magnitude of the reflectivity, the (10), 
(01) and (11) reflections, which peak at the same sample 
rotation angle, are the most likely to be affected. Follow- 
ing Moon and ShuU)^ and taking into account the simpli- 
fication presented by the SANS geometry, one finds that 



to leading order the corrections due to multiple scatter- 
ing are Pw/Pq = Rw + Rio Rii and Ph/Pq = Rh + RIq. 
Here Pio and Pn are the measured scattered power, Pq 
is the measured power of the transmitted incident beam, 
and i?io and Rn are the intrinsic VL reflectivities. Tak- 
ing the peak intensity of the (10) reflection at 0.5 T from 
Fig.[l](w 9 X 10*5 cts./std. mon.) yields Pw/Pq = 0.3%. 
Using the measured Pn/Pw = (0.6)^ w 0.4 (see Fig. [5] 
below) and the above expressions, we find that the differ- 
ence between the normalized scattering and the reflectiv- 
ity are less than 1% for both the (10) and (11) reflections 
and thus insignificant compared to the typical error of 
20% with which the scattering powers can be measured. 

The reflectivity is proportional to the square modulus 
of the VL form factor F{qhk), which is the Fourier trans- 
form at wave vector q^k of the two-dimensional magnetic 
flux modulation of the VL. The reflectivity and the form 
factor for a given reflection is related by 



(1) 



where 7 = 1.91 is the neutron gyromagnetic ratio, t is 
the sample thickness, and qhk is the magnitude of the 
scattering vector;^ Since the vortex spacing, and conse- 
quently q, depends on the magnetic field {q cx \/H), this 
allows |i^((i')| to be determined continuously over wide 
range of scattering vectors using SANS. In the following 
Section we will focus on the field dependence of the first- 
order VL form factor, \F{qio)\^, and in Section UlIDI on 
the higher order form factors. \F{qhk)\'^- 



C. Field dependence of |F(gio)| 

Using the integrated intensity obtained from rocking 
curves, such as the ones shown in Fig. [21 and utilizing 
Eq. ([T]) , one obtains the field dependence of the VL form 
factor of the (lO)-reflections shown in Fig. [3l Within the 
experimental error, field cooling and zero field cooling 
produce identical results; this is indicative of very low 
pining in the sample which is due in part to the post- 
growth annealing.— 

Several models exist for the form factor field depen- 
dence. By far the simplest model is based on the London 
model, extended by a Gaussian cut-off to take into ac- 
count the finite extent of the vortex coresi^i^^ 



Fiq) 



B 



1 + {xqy 



(2) 



Here A and ^ are respectively the penetration depth and 
coherence length, and the constant c is typically taken 
to be between 1/4 and 2^ As shown by the solid line in 
Fig. [21 the measured form factor is well fitted by this 
model which corresponds to a simple exponential de- 
crease with increasing field. Since for all the fields applied 
(Aq)^ 3> 1 the prefactor in Eq. ^ reduces to (j)o/{2nX)'^. 
This is in agreement with our earlier results^ but here 
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FIG. 3: (Color online) Field dependence of the VL (10) form 
factor for both the field cooled (FC) and zero field cooled 
(ZFC) case. The fitted values of the penetration depth and 
coherence lengths are A = 90.7 nm and = 8.22 nm for the 
London model (x^ = 0.14),^^'^ A = 61.9 nm and ^ = 12.7 nm 
for the Clem model (x^ = 0.62), and A = 104.1 nm and 
^ = 7.19 nm for the Hao model (x^ = 1.20).-^'^ 



FIG. 4: (Color online) VL form factor divided by the applied 
field versus scattering vector q for all measured reflections. 
Curves through \F{qio)\ and |-F(gii)| are fits to the London 
model (Eq. ^ and Fig.E]). 



D. Higher order form factors, \F{qhk)\ 



extended to significantly higher fields. Using c = 1/2 
the exponential fit to the form factor yields A = 90.7 nm 
and ^ = 8.22 nm. This value for the coherence length is 
in excellent agreement with previous results. The pene- 
tration depth is about 15% shorter,— which is consistent 
with an improvement of the sample quality by annealing. 
The value of ^ is also in reasonable agreement with the 
estimate based on the upper critical field, ^c2 = 6.7 nm. 

A more rigorous model for the form factor field de- 
pendence was obtained by Clem by including an effective 
core radius and solving the Ginzburg-Landau model.— 
This was later extended by Hao and Clem to include the 
suppression of the bulk order parameter due to vortex 
overlapj '^^'^^ Fits to both of these models are shown in 
Fig. [3l The most noticeable difference between the two 
models is the s-shaped form factor of the Hao model, and 
the significant downturn at higher fields due to the prox- 
imity to the upper critical field Hc2- What is also clear 
from Fig. [3] is that the form factor is somewhat better 
described by the London model (x^ = 0.14) compared to 
either the Clem (x^ = 0.62) or Hao (x^ = 1.20) models, 
and furthermore that the the Clem model returns unre- 
alistic values for A and ^. In this regard, it should also 
be pointed out that the Hao model has in general been 
shown to be a poor approximation to exact, numerical 
solutions of the Ginzburg-Landau modeli^i^ The field 
dependence of the VL form factor in LuNi2B2C, which 
is commonly considered to be a relatively simple super- 
conductor, emphasize the point that any analysis of bulk 
measurements based on a particular theoretical model for 
the VL must be done with the utmost care. Finally, it is 
clear that measurements of F{qi()) alone provides limited 
insight into the VL field distribution. 



The VL form factors for all measured reflections and 
fields are summarized in Fig. [D In both the London 
and Clem model the only field dependence of \F{qhk)\/ H 
comes through the magnitude of the scattering vector, 
Qhk, and therefore the form factors would be expected 
to collapse onto a single curve. This is not observed for 
LuNi2B2C; instead, the form factor follows a different ex- 
ponential field dependence, as shown for \F{qio)\/ H and 
\F{qii)\/ H . Likewise the data does not agree with the 
Hao model, which predicts that for a given q the form fac- 
tor should increase with decreasing field (larger indices h 
and k) and converge towards the value given by the Clem 
model and observed in niobium.— Rather we observe 
that while \F{qii)\/H does indeed he above \F{qio)\/H, 
other higher order form factors fall in-between these two 
limiting curves. 

The deviation from the theoretical predictions is also 
evident if one considers the field dependence of the form 
factor ratio \F{qii)/ F{qiQ)\ shown in Fig. [S] At a field 
of 0.5 T, all the models predict a value of ~ 0.3, decreas- 
ing monotonically by a factor between two and four as 
the field in increased to 4 T. It should however be noted 
that numerical calculations for a square VL in an s-wave 
superconductor predict a largely field independent value 
of \F{qii) / F{qio)\^ which is in agreement with the ex- 
perimental results presented here. 

The failure of the theoretical models to describe the 
measured form factors is not surprising when one keeps 
in mind that they were all derived assuming the screen- 
ing current plane to be isotropic. It is well known that 
LuNi2B2C, as well as the other members of the rare earth 
nickelborocarbide superconductors, posseses a significant 
in-plane anisotropy^'ii^'i^^'i^^ii^i^^ The four- 
fold, in-plane anisotropy manifests itself most strongly in 
the experimental data in Fig. 31 when comparing the [110] 
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FIG. 5: (Color online) Field dependence of VL form factor 
ratios F(qii)/F(gio) (top) and F{q2o)/ F{qio) (bottom). Note 
different scales. The black line is the average and the red line 
is a fit to the data. 



depending on which sign scheme is used. Given that 
the measurements were performed at a temperature and 
field much below both Tc and Hc2, one would expect the 
London sign scheme to apply. This is supported by com- 
paring the field dependence of \F{q2o)/ F{qio)\, shown in 
Fig- El to the numerical work (Fig. 10) by M. Ichioka et 
al. [2|- This shows how the ratio is expected to decrease 
with increasing field, with the sign change separating the 
London from the Ginzburg-Landau regime occuring at 
H K, 1/2 Hc2 for their choice of model parameters. In 
comparison, the ratio \F{q2o)/ F{qio)\ in Fig. [5] extrapo- 
lates to zero at H = l/3i?c2- Further support for the 
choice of all positive form factor signs is obtained by cal- 
culating the magnetic field distribution function shown 
in Fig.[6|c), and comparing the results to muon spin ro- 
tation (/iSR) experiments.— Finally, using the field distri- 
bution (Sneak - -Binin = 1-79 mT) and Eqs. (12) and (13) 
from ref. [8| yields an estimate for the penetration depth 
A = 88.6 nm, in excellent agreement with our fit to the 
London model. 



and [100] crystalline directions, corresponding to the (10) 
and (11) VL-reflcctions. 

IV. DISCUSSION 

While it would be straightforward to incorporate an 
in-plane penetration depth anisotropy into the models 
discussed above, our emphasis here will be on a model- 
free determination of the VL field modulation. 



A. Real-space field reconstruction 

With the VL form factor being simply the Fourier 
transform of the magnetic field modulation, the real 
space field distribution can be obtained from the mea- 
sured form factors by 

B{r)=Y,F{Qhk)e"^'^''-r (3) 

hk 

In the case of the VL the so-called phase problem, aris- 
ing from the fact that only the magnitude \F{qhk)\ is 
measured, is greatly simplified. As the magnetic field 
variation around any vortex exhibits inversion symmetry 
(_B(— r) = B{r)), the form factor must be real and the 
phase problem thus reduces to a sign problem. Within 
the London model the sign on all form factors is expected 
to the be the same, which will be chosen as positive 
and corresponds to having the vortex at the center of 
the unit cell. In contrast, the Ginzburg-Landau model 
predicts both positive and negative signs for the form 
factor, determined by the indices h and k according to 
_(^_l^h^+k^-+hk4i js^Q^g ^Yiat F{qoo) is simply the applied 
magnetic field ^oH. As shown in Fig.[n)Ia,b) the field re- 
construction obtained from Eq. ^ differs significantly 



B. Basal plane anisotropy 

As stated earlier LuNi2B2C possesses a substan- 
tial basal plane anisotropy arising both from the 
Fermi surfac o^'^°i^^'^^ as well as the superconducting 
pairingJ^iiiii^a^iiiaaiiai^a Thermal conductivity mea- 
surements indicates a gap minima, or possibly even 
nodes, along (100),— leading to theoretical speculations 
of a (s -|- g) pairing symmetryii^ Still the dominating 
anisotropy appears to be due to the Fermi surface, as 
indicated by the square VL configuration and discussed 
in Section [TTlSI 

A measure of the in-plane anisotropy can be obtained 
from the SANS results by simply calculating the cur- 
rent flow around the vortices from the field reconstruc- 
tion in Fig.[ni[a) using ^qJ = V x B. Fig.[7]shows | J(r)| 
along the VL nearest neighbor direction as well as the 
VL unit cell diagonal. In cases like this, one frequently 
uses an operational definition of the coherence length (^j) 
as the distance from the vortex center to the maximum 
currenti^i^ From Fig. [7] it is clear that differs for the 
two directions shown. The inset to Fig. [7] shows £,j in 
the basal (screening current) plane. It is striking that the 
minimum is observed along the nodal (lOO)-directions, 
where one would naively expect the weakest pairing and 
hence the largest ^j. It is important to emphasize that 
the field reconstruction is robust in the sense that the 
anisotropy of docs not change (although the absolute 
numerical values do) if one varies the measured form fac- 
tors even well beyond the typical 10% experimental error 
of their values. Rather, our result follows directly from 
the unusually large ratio \F{qii)/ F{qio)\ in LuNi2B2C. 
Finally we also note that the variation cannot be ex- 
plained by a simple "squeezing" effect)^ since the value 
along the nearest neighbor direction ([110]) is larger than 
that along the VL diagonal ([100]). 
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FIG. 6: (Color online) Real space magnetic field reconstruction from the measured VL form factors at 0.5 T and 2 K, using 
the London (a) and Ginzburg-Landau (b) sign schemes. In both cases an equidistant contour spacing of 5 mT was used, 
with the lowest contour at respectively 493 mT (a) and 486 mT (b). In both cases the amplitude of the field modulation is 
48.9 mT or roughly 10% of the average (applied) magnetic field. Four VL unit cells are shown, with a vortex spacing of 64.3 nm 
corresponding to an applied field of fioH = 0.5 T. As discussed in the text the true field reconstruction is with all form factors 
positive (a). Note that the field reconstruction images are rotated 45° with respect to Fig. [1] so that the {110}-directions are 
now horizontal/vertical. Panel (c) shows the magnetic field distribution function for the two different field reconstructions. 



We speculate that our result may share a common 
origin with recent scanning tunneling microscopy mea- 
surements on iso-structural YNi2B2C. These measure- 
ments showed a fourfold-symmetric star shaped vortex 
core extending in the (lOO)-directions at zero energy but 
splitting into 4 peaks and effectively rotating the vortex 
shape by 45° at higher energies within the superconduct- 
ing gap 4^ 

It is interesting to compare the field reconstruction 
and current profiles in Figs.lHIa) and [7] to the numerical 



work of Machida and Ichioka et ai, who have performed 
calculations for a number of VL configurations, pairing 
symmetries and Fermi surface anisotropiesi^^ From this 
it is clear that the best agreement is achieved with an 
anisotropic superconducting gap (d-wave or anisotropic 
s-wave) combined with a Fermi surface anisotropy. Fur- 
ther numerical work to optimize the agreement between 
the experimental and calculated results should provide 
valuable input to calculations of VL configuration, such 
as the one in Ref. [Sll, and result in a more realistic VL 
phase diagram for LuNi2B2C. 



XlO^^ 




f (nm) 

FIG. 7; (Color online) Supercurrent density as a function of 
distance from the vortex center, along the VL nearest neigh- 
bor vortex direction (crystalline [110] axis) and the VL diag- 
onal (crystalline [100] axis). The inset shows the value of 
(distance of maximum current) in the basal plane. To visually 
emphasize the four-fold anisotropy, a circle with radius ^j""' 
is also shown (dashed line). 



V. SUMMARY 

To summarize we have performed comprehensive 
SANS measurements of the VL in LuNi2B2C, thus end- 
ing the common but unsatisfactory practice of discarding 
all but the (1,0) reflection. The measurements confirmed 
the existence of a square VL up to 75% of Hc2- The 
first-order VL form factor, |F(gio)|, was found to de- 
crease exponentially with increasing magnetic field, in 
agreement with the generalized London model but not 
with the supposedly more realistic models for the VL field 
distribution. Measurements of higher order form factors, 
F{qiik); and the real-space reconstruction of the VL field 
modulation, provide a qualitative measure of the in-plane 
anisotropy. This will therefore serve as important input 
to future theoretical work. 

Similar measurement and analysis should be performed 
on other members of the nickelborocarbides; indeed on 
any other superconductor where enough higher order re- 
flections are measurable. In this regard, the results pre- 
sented here will serve as a reference for future work. 
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